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1. Introduction 



1.1 Motivation 

The absence of tree-level flavor-changing neutral currents in the SM means its contributions 
to neutral meson mixing begin at the loop level. The B Q S mixing amplitude is further sup- 
pressed by a combination of the GIM mechanism and the CKM matrix element |Vr s | 2 ~ 2 x 10 
(IV^I 2 ~ 6 x 10~ 5 ), opening the door for observable new physics effects 

A number of recent analyses have revealed hints of new physics in B° s and B° mixing. A 
(2 — 3)<7 tension in the Unitarity Triangle can be explained by new physics in AM S /AM^ The 
UTfit collaboration performed a simultaneous analysis of multiple experimental results, finding 
~ 3a evidence of new physics in B® mixing [||]. D0's recent report of an anomalous like-sign 
dimuon charge asymmetry in semileptonic b-hadron decay deviates from the SM prediction by 
~ 3a, providing additional evidence of new physics in B® mixing [|], |p. 

Disagreement between SM mixing predictions and experiment motivates increased precision 
in the calculation of SM hadronic mixing matrix elements. These calculations have been performed 
on the lattice, with 2+1 dynamical sea quarks, to a (3 — 4)% precision [Q, [7|], but further refine- 
ments are needed to sort out tension with experiment. To permit mixing predictions from new 
physics models, a knowledge of the possible BSM hadronic mixing matrix elements is required. 
Pioneering work [|| [J, utilizing the quenched approximation and static limit of HQET, resulted in 
a quoted ~ 10% precision. We aim to improve upon this. 



1.2 The Role of Lattice QCD 

Whether considering SM or beyond, the disparate scales of hadronization, ^(500 MeV), and 
the underlying flavor-changing physics 1 results in a factorization of the physics at the two scales. 
For example, the SM expression for the oscillation frequency of the B^ <-> 5° transition [10], 

(AM q ) SM = (^|^) ibOO \V tb V t * q \ 2 (B q \O(n)\Bj), (1.1) 

separates the low energy physics of hadronization, characterized by the mixing hadronic matrix 
element from the details of the flavor-changing interactions of the SM. The SM plays a 

role in (B^\0\B 1 ^) by restricting the structure of the effective four-quark operator O. For a generic 
underlying theory we can write 

AM q = Y t c i (n)(B° q \O i {n)\Ej), (1.2) 

i 

where the short distance Wilson coefficients, c,-, contain the details of the underlying theory and are 
generally perturbatively calculable. The ; are all effective, four-quark, AB = 2 (B is the bottom 
quark flavor quantum number) operators allowed by the theory. 

Lattice QCD permits calculation of the non-perturbative, purely hadronic quantities (B^\Oi\B^). 



'in the SM this is the electroweak scale, i€?( 100 GeV), and in new physics models is typically larger. 



2 



and B Mixing in the Standard Model and Beyond: A Progress Report C. Bouchard 

2. Calculation 
2.1 Generating Data 

The space of possible Lorentz-invariant, color singlet, four-quark mixing operators is spanned 



by a basis of mixing operators, referred to for historical reasons as the SUSY basis [11], 

Ox = (b a y^P L q a ) (b P Y ^P L qP), 2 = (b U P L q a ) (f P L qP), 

3 = (b a P L qP) (b P P L q a ), 4 = (b a P L q a ) (f P R qP), and (2.1) 

5 = (b?P L qP)(b P P R q a ), 

listed here with greek color and suppressed spin indices. Of 20 potential mixing operators, 12 can 
be eliminated by Fierz transformation and three by parity symmetry of QCD. We're studying the 
possibility of using these extra operators to effectively increase statistics. For each On we write 
the mixing three-point correlation function as a time-ordered VEV of interpolating and mixing 
operators (a similar, if simpler, process is used to construct two-point correlation functions): 

cTihA) = L(r{(?75%, f2 (Oiv)o ! o(?r5%, fl }), (2-2) 

Xl,X 2 

4s (0 = %{T{{qy 5 by x MY5by Qfl })- (2-3) 

X 

Wick contraction yields products of quark propagators with time ordering ensured by imposing 
?2 > > t\ and t > 0. Heavy (light) quark propagators are obtained by inverting the Fermilab 



(asqtad) action on MILC gauge configurations with 2 + 1 dynamical asqtad sea quarks []12|]. A 
summary of gauge field configurations used (or planned for use) is given in Table [T[ We work in 
the meson rest frame by Fourier transforming the correlation functions and setting p = 0, leaving 



ss a (fm) 


L 3 xT 


P 


mi/m s 


m n L 


N c x N, 


0.12 


24 3 x64 


6.760 


0.1 


3.84 


2099 x 4 


0.12 


20 3 x64 


6.760 


0.14 


3.78 


2110 x 4 


0.12 


20 3 x64 


6.760 


0.2 


4.48 


2259 x 4 


0.12 


20 3 x64 


6.790 


0.4 


6.22 


2052 x 4 


0.09 


40 3 x 96 


7.080 


0.1 


4.21 


1015 x 4 


0.09 


32 3 x96 


7.085 


0.15 


4.11 


984 x 4 


0.09 


28 3 x96 


7.090 


0.2 


4.14 


1931 x 4 


0.09 


28 3 x96 


7.110 


0.4 


5.78 


1996 x 4 


0.06 


64 3 x 144 


7.460 


0.1 


4.27 


826 x N t 


0.06 


56 3 xl44 


7.465 


0.14 


4.39 


800 x N t 


0.06 


48 3 xl44 


7.470 


0.2 


4.49 


668 x N t 


0.06 


48 3 xl44 


7.480 


0.4 


6.33 


668 x N t 


0.045 


64 3 x 192 


7.810 


0.2 


4.56 


860 x N t 



Table 1: MILC ensembles [|12|] planned for use in this study. N c x A^ r is the number of configurations and 
source times (not yet determined if unspecified). This report includes results for the a — 0.12 fm ensembles. 
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Figure 1: (left) Two- and (right) three-point correlation functions in the meson rest frame. The mixing 
operator is placed at the origin in the three-point correlation function. 
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Figure 2: 

the a = 0. 



(left) B s effective mass plot and (right) scaled two-point pseudo-scalar correlation function on 
12 fm, 20 x 64 3 , mi/m s — 0.4 ensemble. Source and sink are lS-smeared. 
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Figure 3: (left) The scaled B s three-point correlation function for mixing operator O4 = (b Lq a )(b Rq^) 
and (right) its relative error on the a = 0.12 fm, 20 x 64 3 , m//m v = 0.4 ensemble. The heavy quark fields 
are lS-smeared at the sinks. 



correlation functions that depend only on time (depicted in Fig. |]). In the three-point correlation 
functions, the heavy quark fields in the mixing operator are improved to remove a discretization 



error [13] and IS -smeared at the sinks to increase ground state overlap. Heavy quark fields in the 
two-point correlation functions are lS-smeared at the source and sink. Figs. ^ and [| show plots of 
two- and three-point correlation function data generated in this way. All data and fit results shown 
are for the B s meson (k = 0.0918 and ma = 0.0349). 

2.2 Fitting 

We build fit functions by decomposing the two- and three-point correlation functions of Eqs. (2.2 
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Figure 4: Simultaneous fits of C 2 ^ and C^ p ' on the a = 0. 12 fm, 20 x 64 3 , m/ jm s — 0.4 ensemble for the B s 
mixing parameter j34, defined in Sect. 3. (left) Fits vs. f m ; n (shown for f max = 24) reveal a common plateau 
for N 3pt = N 2pt = 2 (blue burst), 4 (pink square) and 6 (red dash) fits. Solid lines indicate confidence levels 
for the fits, (right) We plot a representative fit from the plateau (N 3pr — 2 and f^-* = 10) to demonstrate 
stability with respect to f^ax- Fits are blue bursts with error bars and the solid red line is the relative error. 



and |2.3| ) in a basis of energy eigenstates, giving an infinite sum of exponentials 

°° \7 I 2 



n(t+l) ( e -E„t +e -E„{T-t) 



(2.4) 



n,m=0 



2F 

' n \ N\ m } Z l Z m^_y( tl + l) +m ( t2+l ) ( e -En\h\ +e -E n (T-\ti\)\ ( e -E m t 2 +e -E,„(T-t 2 ) 



AE n E n 



where Z n = (n\qy$b). Oscillating opposite parity states, a result of staggered light valence quarks [14], 
and the effect of periodic boundary conditions are accounted for in Eq. (2.4). In practice, we limit 
the time range of data included in the fit ( f min < t < ? max ) and truncate the sums (X^J ~ 1 and ^nm=o ) > 
where N 2 ' 3pt is the number of states used in the fit. 

Fits using data at short times must account for increased excited state contributions by includ- 
ing an adequate number of states. Despite added difficulty, the relatively clean signal in the data at 
short times may make it desirable to include them in the fit. We accomplish this using a Bayesian 



fitting routine [15, 16] and a systematic procedure to select N, t m { a and t max . We are able to achieve 



consistent and stable fits, with a suitable choice of time range, for N 2,3pt = 2,4 and 6. 

From scaled correlation functions we determine time ranges to study. For f m j n we generally 
consider from t m i n = 2 until excited state contributions have significantly decreased. Though in- 
creasing t max utilizes more data, it also introduces an increasing level of noise and can lead to an 
underdetermined covariance matrix. For the two-point correlation function of Fig. ||] we consid- 
ered tj^ = 2,4, 12 and f^ax = 12, 14, ...,24, for the three-point correlation function of Fig. ||, 
tj£ = 2,4, 12 and t^ x = 14, 16, ...,24, and in each case N 2 ^' = 2,4,6. 

We fit for each combination of N 2pt , and tmax an d select a representative fit from the 
common plateau, ensuring stability with respect to our choice of N 2pt , t 2 ^ and f^ax- Then, fixing 
the two-point fit parameters, we repeat the procedure for a simultaneous fit of the two- and three- 
point correlation functions. Fig. ^| demonstrates the stability of the simultaneous fits. 
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2.3 Initial Results 

Table prists preliminary fit results for the B s mixing parameters, /3a?, defined by (B®\On\B®) = 
tf N M Bq ffi (% =2/3,^2 = -5/12,^3 = 1/12,^4 = 1/2, and <T 5 = 1/6). 







mi/m s = 


0.4 


0.2 


0.14 


0.1 


Pi 


X 


{n/afl 2 = 


1.217(14) 


1.196(11) 


1.160(14) 


1.161(14) 


h 


X 


{ri/afl 2 = 


1.509(15) 


1.446(14) 


1.425(16) 


1.448(14) 


ft 


X 


(n/afl 2 = 


1.490(22) 


1.409(25) 


1.318(34) 


1.446(25) 


04 


X 


{r { /afl 2 = 


1.785(14) 


1.731(12) 


1.689(14) 


1.699(13) 


ft 


X 


{r { /afl 2 = 


2.313(26) 


2.255(14) 


2.200(21) 


2.223(13) 



Table 2: Preliminary fit results for on the a — 0.12 fm ensembles. Errors are statistical from the fit. 



3. Outlook 

We are extending B° s fits to other lattice spacings and a range of valence masses, to include B° d . 
We will use these fits in an extrapolation to physical light sea quark mass, the continuum, and light 
valence quark mass (and an interpolation to physical strange quark mass). The continuum [17] and 
staggered lattice [ |l8| ] chiral perturbation theory has been worked out. One-loop perturbative renor- 
malization for the SM mixing operators exists and the BSM operator renormalization is expected 
to be a simple extension of this work. 



Source of Error [%] 


j6i (Lattice 2009) 


Expected 


Reference 


statistical 


2.7 


1.2 


current work 


scale (n ) 


3.0 


1.0 


Oil 


sea & valence quark masses 


0.3 


0.3 




b-quark hopping parameter 


<0.5 


0.1 


[0] 


#PT + light quark discretization 


0.4 


<0.4 


* 


gB'Bn 


0.3 


<0.3 




heavy quark discretization 


2 


~ 1.2 




matching (1-loop perturbation theory) 


~4 


~2.5 


[HI] 


finite volume effects 


<0.5 


<0.1 




Total 


6.1 


~3.2 





Table 3: We estimate an error budget, for J3i, by way of comparison with [|7J. The sources of expected 
improvements are listed in the Reference column. *We anticipate improvement from finer lattice spacings. 

Accurately accounting for errors in the calculation is as important as achieving precision re- 
sults. Our naive statistical errors represent a (40-50)% reduction relative to |[7p. We are generating 
more robust estimates via the bootstrap method. Correlator data exist for a = 0. 12, 0.09 fm and 
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are being generated for a = 0.06, 0.045 fm. Table || quotes statistical and systematic errors for /3i 
from [Q], with expected improvements from the use of a = 0.06 fm data. Analysis of a = 0.045 fm 
data will further improve the error budget. 

Computations for this work were carried out in part on facilities of the USQCD Collaboration, 
which are funded by the Office of Science of the U.S. Department of Energy. This work was 
supported in part by the U.S. Department of Energy under Grants No. DE-FG02-91ER40677 
(A.X.K., C.B., E.D.F.), No. DEFG02-91ER40628 (E.D.F.); the National Science Foundation under 
Grant No. PHY-0555235 (E.D.F.). 



References 



A. Buras, PoS EPS-HEP2009: 024 (2009) [arXiv:0910.1032 [hep-ph]] 



E. Lunghi and A. Soni, Phys. Lett. B 666: 162-165 (2008) |[arXiv:0803.4340 [hep-ph]]| ; J. Laiho, 
E. Lunghi, and R.S. Van de Water, Phys. Rev. D 81: 034503 (2010) 



[arXiv:0910.2928 [hep-ph]]; A. Lenz et al, arXiv: 1008. 1593 [hep-ph] 



M. Bona et al. (UTfit Collaboration), PMC Phys. A3: 6 (2009) [arXiv:0803.0659 [hep-ph]] 



V.M. Abazov et al. (DO Collaboration), Phys. Rev. D 82: 032001 (2010) tarXiv: 1005.2757 [hep-ex]] 



A. Lenz and U. Nierste, JHEP 0706: 072 (2007) [arXiv:hep-ph/0612167] 



E. Gamiz et al. (HPQCD Collaboration), Phys. Rev. D 80: 014503 (2009) 



[arXiv:0902.1815 [hep-lat]]; O. Witzel et al, PoS LAT2009: 243 (2009) 



R. Todd Evans et al, PoS LAT2009: 245 (2009) 



V. Gimenez and J. Reyes, Nucl. Phys. Proc. Suppl. 94: 350-353 (2001) [arXiv:hep-lat/0010048] 



D. Becirevic et al, JHEP 0204 (2002) 025 [arXiv:hep-lat/01 10091] 



A.J. Buras, M. Jamin, and PH. Weisz, Nucl. Phys. B 347: 491 (1990). 

F. Gabbiani, E. Gabrielli, A. Masiero, and L. Silvestrini, Nucl. Phys. B 477: 321-352 (1996) 



[arXiv:hep-ph/9604387] 



A. Bazavov et al. (MILC Collaboration), Rev. Mod. Phys. 82: 1349-1417 (2010) 



[arXiv:0903.3598 [hep-lat]] 



A. El-Khadra, A. Kronfeld, and PB. Mackenzie, Phys. Rev. D 55: 3933 (1997) 



[arXiv:hep-lat/9604004] 



M. Wingate, J. Shigemitsu, CT. Davies, GP Lepage, and H. Trottier, Phys. Rev. D 67: 054505 



(2003) [arXiv:hep-lat/0211014] 



GP Lepage et al, |arXiv:hep-lat/01 10175| (2001). 

GP Lepage's Python-based library for Bayesian least squares fitting of multi-dimensional, nonlinear 
functions of arbitrarily many parameters, built around lsqf it (2008). 



W. Detmold and C.-J. David Lin, Phys. Rev. D 76: 014501 (2007) |[arXiv:hep-lat/0612028] 
J. Laiho and R.S. Van de Water, Collaboration Note (2007). 



A. Bazavov et al (MILC Collaboration), [PoS LAT2009: 079 (2009)1 

R.T. Evans et al. (Fermilab Lattice and MILC Collaborations), in preparation. 

Estimate based on using a = 0.06 fm ensembles. 



7 



